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Can turbulent convective variations drive the Blazhko 
cycle? Dynamical investigation of the Stothers idea 
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ABSTRACT 

The Blazhko-phenomenon, the modulation of the pulsation of RR Lyrae stars remains 
one of the most stubborn unsolved problems of stellar pulsation. The recent idea of 
Stothers proposes that periodic variations in the properties of the convective envelope 
may be behind the amplitude and phase modulation. In this work we approximated 
the mechanism by introducing variations in the convective parameters of the Florida- 
Budapest hydrodynamic code and also by means of amplitude equations. We found 
that the process is only effective for long modulation periods, typically for more than 
hundred days, in agreement with the thermal time scales of the pulsation in RR Lyrae 
stars. Due to the slow response of the pulsation to the structure changes, short pe- 
riod, high amplitude Blazhko-modulation cannot be reproduced with this mechanism 
or would require implausible variations in the convective parameters on short time 
scales. We also found that the modulation of the mixing length results in strong dif- 
ferences between both the luminosity and radius variations and the respective phase 
modulations of the two quantities, suggesting notable differences between the energy 
output of the photosphere and the mechanical variations of the layers. The findings 
suggest that the convective cycle model is not well suited as a standalone mechanism 
behind the Blazhko-cffect. 
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1 INTRODUCTION 

The Blazhko-phenomenon, the amplitude and phase mod- 
ulation of the pulsations of RR Lyrae stars remains one of 
the longest-lived mysteries in astrophysics. Despite of the el- 
egant simplicity it displays at first sight, the Blazhko-effect 
seems to defy any theoretical suggestions. We summarise the 
observational knowledge and compe ting models her e shortly, 
and refer to the excellent review of iKovacsl (J2009J) and ref- 
erences therein for further details. 

Th e Blazhk o -effect was discover e d abo ut a century 
ago by iBlazhkd (|l907l ) and IShaplevi (|l916l ). Since then, 
observations of RR Lyrae stars continued to accumu- 
late, leading eventually to large photometric programs 
as the Konkoly Blazhko Survey led by Johanna Jurcsik 
or The Blazhko Project coord inated by Katrien Kolen- 
berg, using dedicated telescope s (I Jurcsik et al.ll2009bl. ISodorl 
120071 ). multi-site campaigns (IKolenberg et al.l 20061). and 



for RR Lyr see iKolenberg et alj 120061 ) . multiple period 
l|S6dor. Jurcsik fc Szeidll 20101) o r irregularly disappearing 



modulation ( Jurcsik et al.ll2009bl). Longer variations a s the 



inten s ive spectroscopic observ ations ( Chadid fe Chapellierl 
120061 . IKolenberg et all l2010al ). Ground-based efforts un- 
covered numerous interesting details in the Blazhko vari- 
ations of RR Lyrae stars: changing cycle length (e.g. 
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famous "4-year cycle" in RR Lyr l|Detre fc Szeidfl97ot ) were 
reported as well. 

The long and continuous observations of the C0R0T 
and Kepler space telescopes have been providing fur- 
ther details. C0R0T uncove red a modulated RRab w ith 
variable Blazhko-amplitude (jGuggenberger et al.l l201ll ) as 
well as numerous Blazhko-sidelobes up to the 8th or- 
der in the Fourier-spectra of V1127 Aql (|Chadid et al.l 
120101 ). The latter may indicate periodic modulation of 
the amplitude and the phase of the pulsation freq uency 
(|Szeidl fc Jurcsikll2009l . lBenko, Szabo fc Paparoll201ll ). The 
Kepler sample of more than 40 stars revealed many inter- 
esting cases, from hi ghly variable modulation amplitudes 
to p e riod doubling (IKolenberg et al.l 12010b , iBenko et al.l 



l20ld . ISzabo 1 et aljfeoicl . iGuggenberger et ail |2012| ). Addi- 
tional frequencies close to the first and second r adial ovcr- 
tone were also iden tified in both samples |Poretti e t al. 2010, 
IBenko et al.ll20101 ). But from the theoretical side, instead 
of helping to single out a valid theory, all these nuances of 
Blazhko- variables emphasise the shortcomings of the current 
models. 
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1.1 Classical Blazhko- models 



Until the appearance of the Stothers idea IjStothera 12006). 
two models were consi dered as explanations. The mag netic 
oblique rotator model (|Cousenal 19831 . IShibahashill2000T ) pos- 
tulates a strong, ~ 1 kG dipole magnetic field inside the 
star, inclined to the rotational axis. The field distorts the 
radial pulsation, introducing an I — 2 spherical harmonic 
component, and the rotation of the star creates the modu- 
lation pattern, and a quintuplet structure in the frequency 
spectrum. There are two essential drawbacks of the mag- 
netic oblique rotator model: no strong dipo le fields were 
observed unambiguously in RR L yrae stars (IChadid et al.l 
|2004 iKolenberg fc Bagnuldl2009l ). and all the irregularities 
and complexities in the Blazhko-variation (as in the case of 

I| I I 

Jurcsik et al.l 120081 ) contradict with 

the simple geometric, rotation-based explanation. 

The competing model has been the nonradial reso- 
nant rotator which ties the Blazhko modulation to a 1:1 
resonance between the radial and a (preferably 1 = 1 , 
m = 1) nonradial mode (INowakowski fc Dziembowskill200ll . 
iDziembowski fc Cas sisi 1999). Because of the current lack of 
nonradial, nonlinear pulsation models, the resonance model 
was developed using the amplitude equation (AE) formal- 
ism since AEs represent hydrodynamic calculations well. 
The nonradial resonant rotator model, similarly to the mag- 
netic model, proposes a rotation-based mechanism and pre- 
dicts symmetric modulation triplets. Observations are how- 
ever, contradictory: beside the problems with a clockwork- 
like explanation mentioned above, recent, detailed studies 
show that side-peaks are usually neith er symmetric nor lim- 
ited to triplet components ( see e.g. I Jurcsik et al. | | 2005bl . 
iHurtaet al.ll200g| , IS6dorll20ld andlChadid et al.ll2010l). ~ 

The most recent proposal (jBuchler fc Kollathll2011l ) was 
also based on AEs, incorporating the high-order radial res- 
onance th at was found to drive the per iod-doubling phe- 
nomenon IjKollath. Molnar fc Szab6ll201ll ). The calculations 
revealed periodically and irregularly oscillating amplitudes 
that resemble the Blazhko-modulation quite well. These re- 
sults suggest a more intimate relation between radial mode 
resonances and the Blazhko-effect. 



1.2 The idea of Stothers 



The proposed mechanism of lStothersi 1)20061 ) . although lacks 
detailed elaboration, seems to face a number of problems. 
The idea goes as follows: a turbulent magnetic field builds up 
in the outer layers of the star, stalling the convection. Then 
the field is destructed somehow, by ohmic decay or convec- 
tive shredding for example. The periodic changes in the tur- 
bulent convective properties of the star then — in principle — 
are sufficient to change the amplitude and period of the pul- 
sation, creating the observed Blazhko-effect. 

Unlike models that rely on stellar rotation, such mecha- 
nism with stochastic components would allow small changes 
and irregular variations in the modulation. Similarly, mag- 
netic cycles would explain the variations over longer time 
scales, like the 4-year cycle of RR Lyr. A strong but tur- 
bulent magnetic field could explain the lack of detection as 
current observations are well suited only to d ipole fields. But 
there is no mention in the original article |Stothera 120061 ) 
of the required field strength and configuration or how the 



pulsation itself would interact with the mechanism. iKovaca 
(2009) discusses these shortcomings in more detail. The dy- 
namics of the process are also totally omitted as radiative 
and convective single-mode models are compared to each 
other. It is quite possible that on time-scales of the short- 
est observed modulations (a few weeks), the stellar interior 
simply cannot change efficiently if at all to pr oduce observ- 
able m odulation. This was already implied bv lSmolec et al.l 
(J2011T ) . concluding that reproducing modulations, similar of 
observed in RR Lyr, requires huge modulation in the mix- 
ing length parameter. We are interested in a more general 
question: how does the observed variation depend on the 
modulation periods, i. e. how effective is the response of the 
pulsation to the internal modulation of the properties of the 
convective layer? We will investigate these problems both 
with hydrodynamic calculations (Sections [2] and [3| and am- 
plitude equations (Section [4]). 



2 THE MODEL AND THE APPROACH 

General investigation of the proposed mechanism would re- 
quire full three-dimensional MHD modelling of the pulsat- 
ing stars, including radial pulsation, convection and mag- 
netic field evolution. Due to the current lack of such mod- 
els, we explore the overall dynamics of such modulation 
with the Florida-Budapest code, a one-dimensional, turbu- 
lent convective hydrodynamic code. Treatment of the time- 
dependent turbulent conv ection is based on the method de- 
veloped bv lKuhfufj (|l986l ). 

The code involves eight dimensionless a free parame- 
ters of which seven are independent from each other. They 
are connected to the turbulent convective properties of the 
model and are of order unity. Numerical values of a param- 
eters are not provided by theory however, only the compari- 
son with observations provide guidance. We can exploit this 
freedom to introduce modulations in the convective zone. 

The best-known parameter is the mixing length (a\ or 
Q-ml) , and it is often used to fin e-tune the convective proper- 
ties to match the observations. ISmolec et al.l (J2011J) showed 
that large variations in the mixing length produce modu- 
lation and variations in the period-doubling phenomenon. 
However, we usually follow a different path by setting ot\ — 
1.5 as a constant value, and fine-tune the model with the 
other parameters. In those cases we do not consider the mix- 
ing length as an independent parameter. 

A more detail e d des criptio n of the mode l is gi ven in 
IKollath fc Buchler! (l200ll) and IKollath et ail l|2002l) . and 
we also refer to IKollath et alj (|201ll ). where we used the 
same model to investigate the period doubling in RR 



Lyrae stars. We also app ly the amplitude equation method 
IjBuchler fc Gou pil 1984) to describe the modulation char- 



acteristics. 



2.1 Linear results 



Before running the time-consuming non-linear calculations 
to obtain amplitudes, we carried out a li near analysis. Here 
we su mmarise the results published in iMolnar fc Kollathl 
(120101 ) shortly. 

Phase modulation, the change of pulsation period in 
Blazhko RR Lyrae stars usually does not exceed 1-2 %. 
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Table 1. Global parameters of the model and the a parameters 
of the convective zone (mean values for the modulated ones) . 



Figure 1. Linear fundamental period values versus normalized a 
parameters. The large black dot denotes the reference values of all 
four parameters. Dark and light grey areas show the ±0.1% and 
±0.5% changes to the reference periods respectively. Variations 
of pulsation periods of those orders were found in modulated RR 
Lyrae stars. 



Quasi-continuous, space-based photometry allowed to cal- 
culate the instantaneous periods for a few modulated stars: 
four CoRoT RR Lyrae stars reveal ed variations (AP /P to- 
tal amplitudes) between 0.2-1.3% (|Szabo et al.ll2009h while 
V783 Cyg, o bserved with Kepler, sh owed period change 
of ~ 1.1% iJKolenberg et al.l l2010bh . We estimated pe- 
riod changes from published phase mo dulation plots and 
pulsat ion-modulation frequen cy values i n Molnar fc Kollathl 
2010TI: four stars SS Cnc ll Jurcsik et all l200d). RR Lyr 



Kolenberg et all I2006T) AR Her (|Smith et aiPl999l ) and 



MW Lyr 1 Jurcsik et alj 120081 ) have quite similar values to 
each other and to V783 Cyg, between 0.8 — 1 .2% while for 
RR G em (| Jurcsik et al.ll2006r > and DM Cyg jjurcsik et all 
l2009al ) it is somewhat lower, ~ 0.3%. To achieve such 
changes in the linear periods however, quite large variations 
are required in convective parameters. Not all types of a 
parameters are even suitable: for example the eddy viscos- 
ity (ct v ) has very little effect on linear periods since it does 
not change the equilibrium stellar structure. The best can- 
didates we found were the mixing length itself (ax) and the 
parameters controlling the convective flux (a c ) and the tur- 
bulent source function (a s ), but even those would require 
huge modulation amplitudes (Figure [TJ. 

Changes in pulsation amplitudes are related to the 
growth rates of eigenmodes so we examined those dependen- 
cies as well. We approximate the relation from the amplitude 
equation method. The simplest form of an AE with a single 
mode present is A — k,A — qA 3 , where A is the amplitude, 
k is the linear growth rate of the mode and q is the satura- 
tion term respectively. By considering a limit cycle solution 
(A — 0) and a constant q saturation term, one can estimate 
the magnitudes of amplitude variations by A ~ y/H. By ex- 
amining the various relations between the growth rates and 
the a parameters, we concluded that v ariations in a c , a a 
and a „ may be suitable (see figure 2. in iMolnar fc Kollathl 
(J2010J )). But even so, creation of high-amplitude modulation 
requires huge changes that cannot be justified on physical 
grounds in these convective parameters. 
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3 NONLINEAR MODEL CALCULATIONS 

Our goal was to determine whether changes over the typ- 
ical Blazhko-period time scales in the convective environ- 
ment could create suitable modulation in amplitudes and 
phases. We do not attempt to find or validate any underly- 
ing physical processes behind the convective modulation. In 
that sense the modulation int roduced in the m odel is just as 
ad hoc as in original article of lStothersI (|2006l ): some param- 
eters of the stellar structure are varied, and the response 
is observed. There is a crucial difference though: instead 
of comparing unmodulated models with different convec- 
tive parameters, we perturb a limit-cycle solution with time- 
dependent, sinusoidal modulation and observe the variations 
it creates in the global parameters like radius or luminosity. 
We expect that the response of the convective envelope will 

strongly depend on the period of the mo dulation. 

The model in use is the same as in lMolnar fc Ko llath 

(|2010r) . parameters are listed in Table [T] These values place 
it outside the period d oubling instability region (see figure 5. 
in lKollath et al.l 120111 ). Models that reached the limit cycle 
were iterated with modulated a parameters. We considered 
the effects of the mixing length and some individual param- 
eters as well. One reason behind this approach is that tur- 
bulent parameters do not necessarily change in accordance 
in a three-dimensional, interacting convective environment 
as they do in the one-dimensional, mixing length controlled 
models. We particularly emphasize the eddy viscosity pa- 
rameter (otu) as it is strongly connected to the kinetic en- 
ergy through the viscous dissipation and thus to the radius 
variations of the model. Mean radius variations during the 
Blazhko-cycle can be derived from multi colour photometry 
using the I nverse Photometric Method (|S6dor et al.l |2009| . 
IS6dorll2009r ). 

The modulation periods were 10, 30, 100, 300 and 1000 
days, the modulation amplitude of a parameters was con- 
stant in all cases (±25%). The results confirmed our expec- 
tations: the observed modulations in the radius and other 
global parameters grew with longer periods. The hydro- 
dynamic model results are represented with various dots 
on Figures [2] and [3] Note that both the 10 and 30 day 
long modulation resulted in very low amplitude changes in 
most cases. The 10-day period value is cl ose to the shortes t 
modulation observ ed (SS Cnc: 5 . 31 d ll Jurcsik et al.l 12006). 
RR Gem: 7.23 d II Jurcsik et al.l l2005al )) where amplitude 
changes are indeed relatively small: less than 0. 1 magnitude s 
in V colour and about 0.4% in mean radius (|S6dorl 12009). 
The 30-day period is however in the range of many high- 
amplitude Blazhko- variables, including RR Lyr itself (38 d). 
Notable examples with short modulation periods but large 
amplitude variations include MW Lyr, where P m = 15.6 
d, and the mean ra dius variation is calculated to be 1.1% 
IjJurcsik et al.ll2008l ) and the doubly-modulated star CZ Lac, 
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100 
modulation period (days) 

Figure 2. Response to an internal modulation in radius varia- 
tions as a function of modulation period. Squares, dots and cicles 
represent hydrodynamic calculations whereas red, blue and black 
lines represent the amplitude equation results for modulated a„, 
a c and ce\ parameters, respectively. The dashed lines show the 
differences between models calculated with the two extreme val- 
ues of the a parameters, i.e. maximum values the modulation 
may reach (the B(t) forcing amplitude introduced in Eq. [2]l. 



with Blazhko-periods P m \ ~ 18.6 d, P m i ~ 14.5 d, radius 
variations of 0.2 and 0.7 % and brightness vari ations of 0.3 
magn itudes for both modulation components (|S6dor et al.l 
|2010j). These values are quite similar to the results obtained 
for RZ Lyr, a star with significantly longer m odulation com- 
ponents (P m i w 121 d, P m2 » 30 d. lJurcsik et al.ll2012T ). Yet 
the models show higher-amplitude modulation close to and 
above 100 days only, a distinction that does not apply for 
the Blazhko-stars. We note however that for a direct com- 
parison, modulation amplitudes of the radius shall be deter- 
mined along with the variations of the mean properties. 

Not all convective parameters result in the same modu- 
lation properties. For the radius variations (Figure [2]), both 
the mixing length (ax) and the eddy viscosity (a„) display a 
strong modulation period dependence. The convective flux 
parameter (« c ), although it generates only marginal modu- 
lation, nevertheless shows some offset from the correspond- 
ing amplitude equation results (see Section 2]). Differences 
are more pronounced in the case of luminosity (Figure [3]), 
where the modulation amplitude corresponding to the cv c pa- 
rameter is almost constant for all periods and it is notably 
higher, equivalent to the highest <x v values. This difference 
arises from the various effects of the modulations created in 
the model. Coefficients for the amplitude equations were de- 
termined from model series with constant convective param- 
eters. If the convective parameters, especially a c are modu- 
lated, the layer we perceive as the stellar photosphere also 
varies. In addition, the exact determination of the radius 
and the luminosity in the model is not straightforward and 
could result in differences between models with constant and 
dynamically changing convective parameters. This effect is 
also clearly visible for short periods if the mixing length is 
modulated, but for periods longer than about 150 days, the 
hydrodynamic and AE calculations agree: the modulation 
amplitudes increase with longer modulation periods. 

Another feature that differs from the observations is 
the variation of pulsation maxima and minima. Blazhko 
variables show stronger modulation in maxima, our models 



o o 




10 100 1000 

modulation period [d] 

Figure 3. Same as Figure [2] but for luminosity variations. 
Squares, dots and circles represent hydrodynamic calculations 
whereas red, blue dotted and black dashed lines represent the 
amplitude equation results for modulated av, ct c and a\ param- 
eters, respectively. Note the difference between the hydrodynamic 
calculations and amplitude equation results at shorter periods. 



however create more symmetrical variations or sometimes 
even stronger in the minima when transformed to bolomet- 
ric magnitudes. Such discrepancy could be attributed how- 
ever both to the nature of the internal modulation or to the 
limitations of the description of the stellar atmosphere of 
the model. 



3.1 Amplitude versus phase modulation 

An interesting aspect is the relation between the modulation 
of amplitude and phase in the variations of some global prop- 
erties. Here we introduced sinusoidal modulations driven ei- 
ther by the three parameters and compared the properties of 
radius and luminosity variations. In the case of modulated 
a.v, radius and luminosity behave similarly: the magnitude of 
modulation decreases towards shorter modulation periods. 
In other words, the luminosity modulation mostly reflects 
the mechanical variations of the layers of the model. 

If we modulate the mixing length, the changes in radius 
and luminosity variations during the modulation cycle dis- 
play striking differences. The radius variation is relatively 
simple, with decreasing amplitudes and simple curves. On 
the other hand, the phase of the luminosity variation changes 
more wildly and gets stronger towards short modulation pe- 
riods and a looping of the phase relation curve is also visible. 
The modulation properties of radius and luminosity varia- 
tions are displayed in Figured 

These calculations indicate that the modulation of the 
mixing length results in significant differences between the 
oscillating mechanical system and the energy output of the 
photosphere. The effect is strongest if we modulate the con- 
vective flux parameter (a c ) only, but appears in the mix- 
ing length at shorter modulation periods as well. The most 
striking phenomenon is the looping of the luminosity-phase 
curves at 30 and 100 days modulation periods. This fea- 
ture, if confirmed, could provide an additiona l test in itself 
for the convective cycle hypothesis. DM Cyg l|Jurcsik et al.l 
2009a) shows hints of looping but that feature might arise 
from the large scatter of points. Also, the modulation pe- 
riod is 10.57 days, and the model with 10-day modulation 
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Figure 4. Amplitude and phase modulation observed in the ra- 
dius (above) and luminosity (below) variations at different modu- 
lation periods in the model. The large difference between the two 
properties suggest that the photosphere itself experiences strong 
variations if the mixing length (ce\) is varied. Arrows denote the 
direction of progression. Note the looping of the luminosity-phase 
curve at 100 and 30 days. 



(omitted from Figure [3] for clarity) showed no looping. The 
direction of progression is also affected by the looping: long- 
period modulations progress in a clockwise manner but the 
the large loop dominating the 30-day curve and the entire 
10-day curve is purely counter-clockwise. 

We emphasise however that the current description of 
the stellar atmosphere in the model is not well suited to 
investigate the properties of the photosphere in a dynami- 
cally changing environment. A more elaborate analysis with 
a more detailed atmospheric model will be required before 
the observed phase relations (luminosity, radius, pulsation 
phase) and model results can be compared directly. 

3.2 Period doubling 

The continuous and precise data of the Kepler space tele- 
scope revealed period doubling in some modulate d RR Lyrae 
stars (|Kolenberg et al.|[2oTobl . lSzab6 et al.ll2010T ). triggering 
extensive hydrodynamic modelling that traced the origin of 
the phenomenon back to a 9:2 resonance bet ween the fun- 
dame ntal mode and the 9th radial overtone (|Kollath et al.l 
20 111). Period doubl ing was also reproduced in the models of 
Smolec et al.l (|201ll ). Although these results tend to favour 
the mode resonance and convective-cycle models of the mod- 
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Figure 5. Modulated nonlinear model with period doubling. 
Modulation period is 100 days, variations in the a c and a u are 
±75%. Period doubling occurs around and after the peak ampli- 
tude. The time of minimum a values is indicated with the arrow: 
maximum amplitude and period doubling occurs noticeably later. 
Model parameters are: T eff = 6400 K, L = 60 L , M = 0.55 M Q . 



ulation compared to the oblique rotator scenario, neither 
was able to reproduce the modulation it self (without intro- 
ducin g it artificially). On the other hand. lBuchler &: Kollathl 
l)201ll ) incorporated the 9:2 resonance into amplitude equa- 
tions and found not only period doubling (two-mode fixed 
points in the AE formalism) but modulated solutions as well. 

By modulating models that fall into the PD instability 
region, the amplitude variation and/or appearance and dis- 
appearance of the PD phenomenon during the modulation 
cycle can be compared to the observations. On one hand, 
data from the first 127 days of Kepler (Ql and Q2) indi- 
cate that period doubling is the strongest on the descending 
or ascending branches of the Blazhko-modulation, and vari- 
ations from one B lazhko-cycle to another are also obvious 
jSzabo et al.ll20ld ). On the other hand, period doubling oc- 
curs predominantly during th e Blazhko max i ma in the mod- 
els, as already pointed out bv lSmolec et al.l (|201lT ). and the 
cycles are repetitive, as expected from a repe titive modula- 
tion. The inspection of Floquet stability roots (jKollath et al.l 
l201ll ) also showed that instability is expected towards less 
efficient convection (smaller a parameters, closer to the ra- 
diative limit) where the pulsation amplitudes are higher. 

A nonlinear model with amplitude modulation and pe- 
riod doubling is shown in Figure(5] Parameters are: 0.55 Mq, 
601/0, 64007^, modulated parameters are a c and a v , both 
by ±75% with a period of 100 days. Unusually, period dou- 
bling appears around the peak amplitude and dominates the 
descending branch of the amplitude modulation. The effect 
is very similar if the mixing length itself is modulated. 



4 AMPLITUDE EQUATION CALCULATIONS 

Nonlinear model calculations are time-consuming exercises 
especially if long-period modulation is involved. If the time 
scales are separated enough, amplitude equations can be 
used instead of the hydrocode, excluding the pulsation from 
the in tegration and focusing only on the amplitude varia- 
tions (jBuchler fc Goupillll984l ). Coefficients have to be de- 
termined first from the models but the freedom of changing 
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them is another possible advantage. The effects of different 
growth rates for example can be investigated much faster 
this way. 

As noted in Section 12.11 the simplest amplitude equa- 
tion with only one mode present is A — k,A — qA . The 
coefficients of AEs were determined the following way. Lin- 
ear growth rates of the fundamental mode («(«)) were calcu- 
lated in the linear models with different a parameters. Then 
the corresponding nonlinear models were iterated and satu- 
ration coefficients were determined from constant amplitude, 
limit cycle models (A = Ao), in the form of q(a) — k/Aq. 
We considered variations in three model parameters, a.,,, a c 
and a>\. For all parameters, the equation: 



A(t) = «(a(t)) A - q(a(t)) A 3 



(1) 



had to be solved. By setting n/q = B 2 , the equation can be 
expressed as: 



A ( A 2 



(2) 



where B is the forcing amplitude exerted by the changes of 
the stellar interior: without the non-zero reaction time of the 
pulsation, the observable amplitude would be B(t). Ampli- 
tude changes are however governed by the above equation 
which is a nonlinear filter with phase shift and a frequency 
response defined by k. To create large variations over shorter 
time scales, B has to be similar or even higher than A which 
is physically implausible. These properties are shown by the 
integration of fTJ in the following subsections. 



4.1 Response to an internal modulation 

We calculated AEs with a constant modulation amplitude 
(a„, a c and a\ varied by ±25%) and a wide range of modu- 
lation periods. The observed modulation in the stellar radius 
is plotted in Figure [2] The solid lines show the AE results 
while hydrodynamic calculations are represented with vari- 
ous dots for comparison. Amplitude equations represent hy- 
drodynamic calculations well in general, however some dif- 
ferences arise in our case. As we explained it in Section [3l 
these effects are caused by the inherent differences between 
the calculated luminosity and radius values of constant and 
modulated hydrodynamic models. Even with these restric- 
tions, the results clearly show that short and even medium- 
period modulation (< 100 days) results in low-amplitude 
modulation of the observable quantities. The highest pos- 
sible amplitudes (i.e. the difference between the two ex- 
treme states, dashed lines) are reached only at very long 
modulation periods, over 1000 days. The values represented 
by the dashed lines are similar to the 8R/R value used by 
IStothera (|2011f ). the difference between the radial and con- 
vective models. 

In all three cases the observed modulation amplitudes 
increase linearly from short to long modulation periods be- 
fore approaching the maximum value asymptotically. The 
only difference is in the highest value: modulation of the 
convective flux has the smallest effect on the radius vari- 
ations followed the eddy viscosity and the mixing length. 
The luminosity values are equally affected by the former 
two parameters but both create about the third of the effect 
of mixing length only. Modulating the system with larger 
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Figure 6. Phase difference between the modulation of the a v 
parameter and the radius variations versus the period of the mod- 
ulation. Values approach tt/2 for moderate and fast modulation. 
Inserts show two cases: the modulation period for the upper right 
panel is P = lOOOd while for the lower left panel P = 31d Blue 
dashed lines show the B(t) forcing amplitudes while the solid 
black lines show the A(i) response amplitudes. The phase dif- 
ferences between the maxima of the two variations are indicated 
with the black bars. 



amplitudes allow for larger observed variations on shorter 
time scales as well but in that case even stronger modula- 
tion should occur over long periods. Of course the possibility 
of slow but very strong modulation cannot be ruled out by 
the lack of observations. A much stronger argument against 
such mechanism is the existence of fast, strong modulation 
in some RR Lyrae stars: it would require enormous changes 
in the convective properties of the stellar envelope over only 
a matter of weeks. 



4.2 Phase relations 

The slow response of the pulsation to the internal modula- 
tion also manifests itself as a phase shift between the internal 
modulation cycle and the observable variation of the pulsa- 
tion amplitude. The phase difference is shown in the main 
plot of Figure [6] For modulation periods under about hun- 
dred days the difference approaches 7r/2. Two examples are 
shown in the subplots with modulation periods P — 1000 
and P = 31 days. Solid and dashed lines represent the A(t) 
response (observable) and B(t) forcing amplitudes as defined 
in @. The decrease of amplitudes with shorter modulation 
periods is also very prominent. These results are confirmed 
with the corresponding hydrodynamic calculations as well. 
Figure [S] shows that both the highest-amplitude pulsation 
and period doubling occur after the a parameters reached 
their minimum values. 



4.3 Changing the growth rate of the mode 

The normalised linear growth rate of the fundamental mode 
in RR Lyrae models is typically in the range of one per cent. 
Amplitude equations allow us to change these values eas- 
ily and compare cases with different growth rates. Figure [7] 
shows the effect of multiplying the growth rates by 2 and 10. 
Since Aq = KJq in our simple model, the saturation term was 
multiplied as well to fix the magnitude of the amplitudes. 
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Figure 7. Changing the growth rate of the mode shifts the curve 
along the time axis. Higher growth rates would allow faster mod- 
ulation as well. The red solid line (k) is the same as the red solid 
line («„) on Figure [2] The blue dashed (2k) and black dotted 
(10k) lines show the effect of higher growth rates. 



Varying these parameters essentially scales the time param- 
eter to an arbitrary unit. But these simple exercises further 
confirm that the slow response of the pulsation mode is the 
main drawback for the Stothers-model: the 10k (tenfold in- 
crease in the growth rate) case allows much faster amplitude 
changes: large amplitude variations over 10-20 pulsation pe- 
riods, or a few dozen days. In contrast, iBuchler fc Kollathl 
i|201ll ) showed that resonant coupling with the 9th overtone 
can generate adequate amplitude modulation with reason- 
able mode growth rates. 



5 CONCLUSIONS 

In this paper we i nvestigated some dynamical aspects of the 
idea proposed by IStothersI (|200q ) to explain the Blazhko- 
effect. Although the advocated scenario is very complicated, 
the idea is based on the simple comparison of some prop- 
erties (amplitudes, periods) of radiative and convective RR 
Lyrae models. The properties of the mechanism, involving 
a highly variable turbulent magnetic field coupled with the 
convective properties of the envelope, represent great chal- 
lenges both for observations and detailed stellar models. We 
restricted ourselves to a quite simplified method by planting 
an ad hoc internal modulation into a ID turb ulent convec- 
tive m odel. A similar approach was followed bv lSmolec et al.l 
i|201ll ) to recreate the modulation properties of RR Lyr. 
However, instead of focusing on the best reproduction of 
the modulation in a single star, we carried out a more gen- 
eral analysis, studying the possible modulation amplitudes 
over a broad range of modulation periods. 

As the mean pulsation period also changes during the 
Blazhko-cycle, we first compared the periods of linear mod- 
els with different convective parameters. The two suitable 
parameters to achieve the required variations (^ 1%) were 
the ones controlling the convective flux and the turbu- 
lent source function but even those require huge variations. 
These results already present restrictions for the further, 
non-linear calculations. 

Blazhko RR Lyrae stars come in all flavours, most im- 
portantly in all kinds of different modulation periods and 



amplitudes. The greatest challenges for the mechanism are 
the ones with strong amplitude variations over short time 
scales, under about 40-50 days. We investigated the dy- 
namics of the mechanism by modulating various convective 
parameters (eddy viscosity, strength of the convective flux 
and mixing length) in the model and by comparing the ob- 
served amplitude variations of the global stellar parameters 
over different modulation periods. We found that the mech- 
anism is efficient only for long modulation periods (> 100 
days) but fails to explain the short-period, large-amplitude 
Blazhko stars. Increasing the internal modulation also pos- 
tulates huge amplitude changes for very long modulation 
pe riods. These r e sults are in agreement with the findings 
of ISmolec et al.l i|201ll ). It is worth mentioning that the 
ch anges applied to con vective parameters in this paper or 
bv lSmolec et al.l 11201ll) ar e still less than those in the orig- 
inal paper (|Sto thcrs 2006J) where simply a fully convective 
and a fully radiative model were compared. 

The reason of the ineffectiveness of the mechanism is 
the slow response of the pulsation to the changes in the stel- 
lar envelope. Typical normalised growth rates of the modes 



in RR Lyrae stars are in the range of ~ 10 



10" 



The 



mechanism of Stothers would require modes with an order of 
magnitude faster growth rates to allow efficient amplitude 
changes over a few dozen pulsation cycles or less. The re- 
sults also suggest that although the difference between the 
convective and radiative mod els seem to acc ount for the ob- 
served mean radius changes (|Stothersll201ll ). the process is 
dynamically incapable to reproduce the variations under the 
modulation periods mentioned in the paper {e.g. 15.6 days 
for MW Lyr). The actual difference between the two models 
is reached only when incorporating years-long modulation 
periods. 

The results also indicate that the luminosity and thus 
the properties of the photosphere may exhibit variations 
which are different from the radius that reflects the changes 
in the pulsation energy only. We found that the observed 
luminosity modulation is less dependent on the modulation 
periods but even in the best case, only moderate modulation 
could arise for modulation periods under about 150 days. We 
note however that a more detailed description of the stellar 
atmosphere is needed to draw further conclusions {e.g. phase 
relations) about the properties of the photosphere in a dy- 
namically changing environment. 

All these results suggest that the mechanism proposed 
by Stothers faces at least as serious shortcomings as any 
of the current Blazhko-models. Most prominently, it is not 
suitable to explain strong modulation on short time scales 
especially without involving large amounts of internal vari- 
ations in the convective parameters. Our analysis of course 
neglects a number of processes, e.g. the feedback of pulsation 
to the magnetic field and the dynamical coupling between 
the convective envelope and the variable magnetic field is 
only crudely approximated with the internal modulation. 
Without such a detailed modelling effort, the plausibility of 
the model cannot be rejected. Instead of a standalone mech- 
anism though, it may be incorporated into a more complex 
explanation of the Blazhko-effect along with the recent find- 
ings about period doubling and radial mode resonances. 
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